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Abstract. We present detailed calculations of low-energy spin dynamics in the "ferric wheel" systems 
Na:Fe6 and CsFeg in a magnetic field. We compute by exact diagonalisation the low-energy spectra and 
matrix elements for total-spin and Neel-vector components, and thus the time-dependent correlation func- 
tions of these operators. Comparison of our results with the semiclassical theory of coherent quantum tun- 
nelling of the Neel vector demonstrates the validity of a two-state description for the low-energy dynamics 
of ferric wheels. We discuss the implications of our results for mesoscopic quantum coherent phenomena, 
and for the experimental techniques to observe them, in molecular magnetic rings. 

PACS. 75.10.Jm Quantised spin models - 03.65.Sq Semiclassical theories and applications - 73.40. Gk 
Tunnelling - 75.30.Gw Magnetic anisotropy 



1 Introduction 

Quantum coherent tunnelling of the magnetic moment 
in nanoscopic magnets has recently become the focus of 
strong experimental and theoretical activity M . The ferric 
wheel systems Fejv (Fig. 1) present a particularly promis- 
ing subgroup in which crystals have now been prepared of 
compounds with N = 6, 8, 10 ,12 and 18 magnetic Fe(III) 
ions in ring geometry 

@ISill0>i- These molecules have 
antiferromagnetic (AF) coupling between spins s = 5/2 on 
each iron site, show a ground state with vanishing total 
spin S = at zero field, and because of an effective uni- 
axial magnetic anisotropy admit the possibility of meso- 
scopic quantum phenomena in the form of coherent tun- 
nelling of the Neel vector §,0|ll|,|l|. 

The best characterised molecular rings are Feio J|, 
q|, various realisations of Fe6 which differ in lig- 

and group and central alkali metal ion, and Cs:Fes j§. 
These materials have been studied by a variety of exper- 
imental techniques, includi ng mag netic susceptibility and 
torque magnetometry [§,|l3], jj .jig], specific heat ]l(|, elec- 
tron spin resonance (ESR) [17|, inelastic neutron scatter- 
ing (INS) jjS and spin relaxation in nuclear magnetic res- 
onance (NMR) . All of these studies serve essentially to 
characterise the zero-field spectrum and the dependence 
of the energetic separation of the lowest two levels on the 
applied field, both of which may be encapsulated within 
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Fig. 1. Schematic representation of an Fe6 ring with spins 
aligned along the easy axis (z), whose orientation is normal to 
the ring plane. The held h = g/xsB is applied in the plane of 
the ring. 



a phenomenological Hamiltonian written in terms of only 
the total spin f||l5|]. 

In contrast to the situation in the ferromagnetic (FM) 
molecular clusters Mni 2 and Fes the notion of spin 
quantum tunnelling in the AF ferric wheels has to date 
received little experimental attention. On the theoretical 
side, a semiclassical description of the low-energy dynam- 
ics provides the clear prediction jl^] of coherent tunnelling 
of the staggered moment. Analysis of magnetisation and 
torque measurements in this framework |2Cfl allows one to 
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extract the AF superexchange J and an effective uniax- 
ial anisotropy k z for the ring molecules. k z determines the 
height of the tunnel barrier between the degenerate energy 
minima, and thus the extent to which quantum transitions 
may influence the system response at the lowest temper- 
atures. 

Although technical difficulties certainly arise in cou- 
pling directly to the staggered moment, the ferric wheels 
appear to offer significant advantages over FM molecules 
for the experimental observation of coherent tunnelling. 
In all such molecular crystals there are many possible 
sources of decoherence |20] which act to destroy the co- 
herent nature of the predicted tunnelling processes. These 
we discuss in further detail in Sec. IV. Here we note only 
that any measurement of quantum coherence effects would 
require at minimum that the decoherence time r^ 1 be 
significantly longer than the tunnelling time Z\ -1 . From 
the semiclassical theory, summarised in Sec. 0, A is the 
level separation in the low-energy (two-state) manifold, 
and corresponds to the tunnel splitting. In the FM sys- 
tems f»zl, and mesoscopic spin quantum tunnelling is 
said to be incoherent. While the reasonable assumption 
that decoherence rates, r, are not vastly different in the 
FM and AF systems remains to be proven, the tunnelling 
frequency, A, is some 6-7 orders of magnitude larger in 
the AF ferric wheels, making these very good candidates 
for the observation of coherent tunnelling. We note here 
that none of the above types of experiment offers a means 
to extract a decoherence rate or to distinguish between 
coherent and incoherent tunnelling. 

The decoherence rate can be determined only from dy- 
namical quantities plj ] , whose spectral linewidths provide 
an upper bound on r. In the strictest sense, experimen- 
tal confirmation of the inequality r <C A still does not 
establish the existence of coherent dynamics, an under- 
taking which would require a true time-domain observa- 
tion of an appropriate oscillating quantity. In this study 
we seek to establish that coherent oscillations, including 
coherent tunnelling processes, are indeed present in the 
dynamical properties of ferric wheel systems without de- 
coherence. That this is not a trivial statement is clear both 
phenomenologically from the apparent success of a model 
requiring only the total spin H,[L3j, and microscopically 
from the huge number of states |j2s + 1)^], energy-level 
splittings and possible matrix elements involved in a com- 
plete description. Establishment of the presence and na- 
ture of quantum coherence in these mesoscopic molecular 
systems would provide both an existence proof for coher- 
ent spin tunnelling if the condition r <C A is satisificd, 
and valuable guidance for its experimental observation. 

Thus we present a detailed investigation of the dy- 
namical properties of the ferric wheels Na:Fe6 and Cs:Feg 
by exact diagonalisation (ED). A necessary initial step 
is to identify those dynamical quantities from which in- 
formation on the quantum dynamics of interest may be 
obtained. This task is straightforward for FM molecular 
clusters, but is less immediately evident in ferric wheels, 
where, with the exception of a recent examination of elec- 
tronic and nuclear spin dynamics in a semiclassical frame- 



work 21 , a full analysis of microscopic dynamical prop- 
erties is still lacking for the realistic Hamiltonian. The 
spin dynamics of rings with AF Heisenberg interactions 
has attracted some recent interest (see Refs. [^2,23 and 
references therein). However, previous studies have been 
restricted to systems without anisotropy, which are not 
expected to show mesoscopic quantum phenomena in the 
form of coherent Neel vector tunneling. 

In contrast, we consider here the spin dynamics of a 
full, effective Hamiltonian for the ferric wheels Fe6 and 
Fes- The calculation of dynamical quantities, which we ob- 
tain from total-spin and Neel- vector correlation functions, 
then provides essential new information adding to the un- 
derstanding of the physical properties of ferric wheels. 
Such a study is required not only to establish the exis- 
tence of quantum coherent oscillations in a microscopic 
model for a complex, mesoscopic system, but also to aid 
the extraction of decoherence rates from experimental dy- 
namical quantities. On the theoretical level, dynamical 
studies are qualitatively more difficult than the calcula- 
tion of ground-state properties Q because they require 
that excited states and all corresponding matrix elements 
be taken into account. We note that a very recent, related 
study H] of the systems Fe6 to Fex2, while technically 
advanced, includes the anisotropy only at the level of a 
low-energy effective Hamiltonian and does not consider 
intrinsic dynamical correlations in the presence of a time- 
independent magnetic field. 

The manuscript is organised as follows. In Sec. || we 
present the model Hamiltonian and observables, and pro- 
vide a brief overview of the technique by which they are 
analysed. In Sec. || we present our results for Na:Fe6 and 
CsiFes, both in the time domain and by analysis of matrix 
elements, and illustrate their physical origin by compar- 
ison with semiclassical approaches. In Sec. ^ we discuss 
the microscopic understanding of quantum coherence, de- 
coherence sources and the prospects for experimental ob- 
servation of spin tunnelling in ferric wheels. Sec. |5| contains 
a summary and conclusions. 



2 Model and Method 

We work within the minimal model Hamiltonian for AF 
rings with effective uniaxial anisotropy ]l^,p0|[21[|, 



N 

H = J Si s 



JV 



N 



(1) 



where N = 6 or 8, si — s^+i and h = g/i^B. J is the su- 
perexchange interaction which favours an AF spin config- 
uration on the ring, k z is the effective uniaxial anisotropy, 
of dipolar and single-ion origin, and the final term is the 
Zeeman coupling. In the following all energies and fields 
are scaled to J and h is set to unity. We define the total- 
spin operator 

N 

S = $> (2) 
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and the Neel-vector operator 



1 N 
Ns ^ K ' 



(3) 



For finite anisotropy the total spin is not a good quantum 
number, meaning that k z leads to mixing of different spin 
multiplets, as will be evident from anticrossings in the en- 
ergy spectra as a function of field to be shown below. This 
is a principal reason why numerical calculations are re- 
quired to make further progress in a quantitative analysis 
of Eq. (|l|). The dynamical variables we consider are the 
autocorrelation functions of the total-spin and Neel-vector 
operators, 

S aa (t) = (S a (t)S a (0)) (4) 



and 



Af aa (t) = (n a (t)n Q (0)), 



(5) 



from which one may seek temporal oscillations character- 
istic of coherent tunnelling dynamics. Working with cor- 
relation functions of total spin and Neel vector is advan- 
tageous because it allows one to retain some of the spatial 
symmetries. This simplifies the computation of matrix el- 
ements of S and n required in addition to the energy spec- 
tra for a full understanding of the dynamics. 

The dynamical response of a single spin pj| may also 
be considered directly. However, we state that all relevant 
dynamical properties are encoded in the two correlation 
functions above [Eqs. (|^,|^)], and show that the single-spin 
quantities can be deduced from these as follows. We denote 
by e lp the eigenvalue of the one-site translation operator 
on the ring, and presume that all low-energy states are 
contained in the sectors p = and p = n, a fact we will 
verify below. If the ground state, |0), is in the sector p — 0, 
then matrix elements (i|S a |0) are finite only for states \i) 
in the sector p — and elements (i|n Q |0) are finite only for 
states \i) in the sector p = tt. Because all states considered 
are invariant under translations by two lattice sites, S may 
be substituted by |iV(si + s 2 ) and n by t^(si — s 2 ). It 
follows that single-spin correlation functions are given at 
low frequencies by 



1 



(s la (t)s la (0)) ~ sW na (() + j^S aa (t), 



(0) 



because the cross-correlation functions (n a (t)S a (0)} van- 
ish due to the opposing symmetries of n and S under 
one-site translation. 

Taking as a guide the semiclassical treatment of Ref. 
||l2| , a tunnelling scenario is applicable if the ground and 
first excited states, |0) and |1), are energetically well sep- 
arated from all other states and form a (weakly) tunnel- 
split doublet with splitting A = E\ — Eq. Because the 
states |0) and |1) have opposite behaviour under trans- 
lation, as defined above, the total-spin matrix element 
vanishes in this manifold, (1|S|0) = 0, and S zz {t) has 
no coherent oscillations with period characteristic of the 
tunnelling time By contrast, the dynamical prop- 

erties of the Neel vector are dominated by tunnelling in 
the lowest manifold, |(l|n z |0)| ~ 1, so that Af zz (t) should 



exhibit coherent oscillations with period 2ir/A. The quan- 
tities S and M are measured through the susceptibili- 
ties xs and xn, which are directly related to Eqs. (|]]|) 
by the fluctuation-dissipation theorem. Xs ^ s accessible 
in ESR or alternating-current (AC) susceptibility mea- 
surements, but, because Xs( w ^ A) = shows no re- 
sponse at the splitting frequency A, does not contain any 
information on tunnelling dynamics pl| , p5| , by which is 
meant processes in the low-energy sector. is the quan- 
tity which, in accordance with semiclassical theory [ [i"2[ 
BjJ , should show oscillatory behaviour due to coherent 
Neel vector tunnelling. In the weak tunnelling regime, 
%^(w ~ A) ~ tt8(uj — A) tanh(/3Z\/2) has a delta-function 
peak at A plf , but its experimental observation is not 
straightforward. 

A microscopic analysis of Eq. ([!]) involves many en- 
ergy levels coupled by potentially large matrix elements of 
the operators S and n arising from the spin interactions. 
Qualitatively, the system in a transverse magnetic field 
(Fig. 1) exhibits two degenerate, classical spin configura- 
tions (obtained by reflection in the ring plane), between 
which the semiclassical approach predicts a tunnelling sce- 
nario. However, because of the approximations involved in 
this description, it is not clear that a fully quantum me- 
chanical treatment of the ferric wheel systems would con- 
firm the presence of two low-lying levels sufficiently well 
separated from all others, or that Af aa would exhibit co- 
herent oscillations dominated by a single frequency for any 
choice of a or of the applied field. We use the exact corre- 
lation functions, meaning the spectra and matrix elements 
required in their calculation, to resolve this issue. 

Magnetisation curves of Heisenberg ring s (Eq. (0) with 
k z = 0) were computed some time ago using Lanczos ED 
[26 1, while recent computations of dynamical properties 
[23 1 remain very similar in scope. The situation becomes 
more complicated in the presence of a non-zero single- 
ion anisotropy k z ^ 0, and of a magnetic field applied 
at an arbitrary angle to the z-axis. First, in this case 
the eigenvectors have a non-trivial dependence on the ap- 
plied field. Second, when S z is no longer a good quan- 
tum number, the loss of the associated symmetries in spin 
space causes a substantial increase in the dimensions of the 
Hilbert spaces. With modern computers it is possible to 
compute ground-state properties by the Lanczos method 
without exploiting spatial symmetries, at least in the case 
of Fe6 [p0[ . However, as indicated above, dynamical stud- 
ies require in addtion many excited states and the corre- 
sponding eigenvectors, although some spatial symmetries 
remain which simplify this task. The Lanczos method can 
also be used for this purpose J27j], but becomes substan- 
tially more involved. 

The complete Hilbert spaces for Fe6 and Fes rings con- 
tain 6 6 = 46656 and 6 8 = 1679616 states, respectively. For 
the general case of Eq. (0), there are no symmetries in spin 
space |2£| , but the spatial symmetries of one-site transla- 
tion and reflection at a given site are present. The ground 
state in a magnetic field is always located either in the sec- 
tor with p = or with p = 7r, and has positive parity with 
respect to reflection at the selected site. These two sec- 
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tors also contain the lowest excitations. The dimensions 
of these sectors are 4291 (107331) for p = and 4145 
(106680) for p = tt in the case of Fe6 (Fes). Since the Neel 
vector connects these two subspaces, we work in the sum 
of the two spaces, which has dimension 8436 (214011), for 
the computation of matrix elements and correlation func- 
tions. 

There are several viable methods for the computation 
of a large number of extremal eigenvalues and -vectors. We 
have used a combination of simultaneous vector iteration 
of a large number of vectors with explicit diagonalisation 
of the Hamiltonian H in the subspace spanned by these 
vectors during the iteration (29). For Fe6 we have com- 
puted the lowest 350 eigenvectors and eigenvalues in the 
symmetry subspace described above. The error caused by 
this truncation may be estimated in two ways. One is to 
examine how the results change on varying the number 
of states retained, and the other is an exact evaluation at 
t = of the correlation functions (jij) and (^) for compari- 
son with the results obtained using the truncated spectral 
representation. On the basis of both methods we estimate 
that the truncation leads to errors in the temporal correla- 
tion functions for Feg on the order of 10 -6 of the peak val- 
ues, meaning that the truncation error is undetectable in 
any of the figures to be shown below for Feg, while numer- 
ical errors for the individual eigenstates are considerably 
smaller still. For Fes we have retained the lowest 50 eigen- 
vectors in the combined subspace, and have an estimated 
truncation error on the order of 10~ 4 of the peak values 
in the temporal correlation functions. As noted above, all 
of our results are obtained for zero temperature using the 
minimal Hamiltonian ([!]), and thus include no sources of 
decoherence. 

As in Ref. |m| we will focus on magnetic fields ap- 
plied in the plane of the ring (Fig. 1), the geometry which 
retains the highest tunnel barrier, with maximal locali- 
sation of the Neel vector, at given field. We have also 
considered field angles out of the plane of the ring, and 
confirmed that they yield qualitatively similar results at 
strong fields. This is to be expected because the physi- 
cal situation remains one in which the spins are largely 
confined to a planar motion between two potential min- 
ima, and also because in our numerical approach the z- 
inversion symmetry of the transverse-field case is not es- 
sential. Here we calculate dynamical quantities [Eqs. (^j,||)] 
exhibiting oscillations in the time domain, and discuss the 
information they contain concerning mesoscopic quantum 
coherent phenomena. 

3 Results 

3.1 Na : Fe 6 

We begin with a discussion of the energy spectrum for 
the physical system Na:Fe6, for which k z /J — 0.0136 in 
Eq. (|l]) |20). As the applied field is increased, level cross- 
ings occur at critical fields B cn between ground states with 
increasing total spin and alternating quantum numbers 
p = or p = n. This leads to a magnetisation curve with 



T 
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Fig. 2. Lowest energy level spacings Ai = Ei — Eq for an 
Na:Feg ring described by the minimal Hamiltonian (jl|), for ap- 
plied fields gfiBBx = 3.1J (a) and g^BB x = 3.5J (b). The 
thickness of the vertical lines represents the magnitude of the 
matrix elements of n z connecting each level pair. 

an almost regular staircase of plateaux p3] . Guided by the 
semiclassical prescription that quantum tunnelling is best 
defined in intermediate fields |l2j| , we consider magnetic 
fields beyond the lowest magnetisation plateau (B > J3 c i). 

Fig. 2 shows the lowest energy levels for two fields cho- 
sen near the centre of a magnetisation plateau (a) and very 
close to a level crossing (b) (see Figs. 5 and 6 below). Near 
the level crossing, there are indeed two nearly degenerate 
levels lying well below any of the others, but we stress that 
this alone is not sufficient to guarantee single- frequency os- 
cillations corresponding to the energy difference A in any 
observable, and thus to justify a two- level tunnelling sce- 
nario. This point is represented schematically in Fig. 2 by 
the lines connecting the levels, the thickness of which cor- 
responds to the magnitude of the matrix element |(i|n 2 |j)|, 
where i and j denote the energy levels. For fields corre- 
sponding to the centre of a plateau, one observes by con- 
trast that there is no clear two-level manifold, but that 
for certain operators, such as n z as shown, the matrix el- 
ement between the lowest pair of levels is dominant. This 
situation, |(l|ri z |0)| ~ 1, can be taken to express the re- 
quirement for a two-level description to be adequate for 
mesoscopic tunnelling of the Neel vector. 



3.1.1 Time Domain 

There are in principle two ways of testing the coherent 
low-energy dynamics of a two-level system at low tem- 
peratures. The first would be to prepare the system in a 
non-eigenstate of the Hamiltonian, \ip) = (|0) + |l))/\/2, 
and then to observe coherent oscillations of the quan- 
tity of interest (here n z ) in the time domain. The second 
is to measure ground-state correlation functions such as 
(0\n z (t)n z (0)\0). For an idealised tunnelling scenario in 
which n z connects only |0) and |1), these two quantities 
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Fig. 3. The Neel vector correlation function N zz (i) for Na:Fes 
in magnetic fields (a) g\j,sB x = 3.1J and (b) g\isB x = 3.5J. 





B x [J/g^ B ] 



Fig. 5. Evolution with magnetic field of (a) the energy-level 
splittings Ai and (b) the matrix elements (i|n z |0) for the low- 
energy sector (i = 1,2) in Na:Feg. ED results (symbols) are 
compared with the RR approximation (lines). 




t[1/J] 

Fig. 4. The total-spin correlation function S zz {t) for Na:Fe6 
in magnetic fields (a) g\j,sB x — 3.1J and (b) g\isB x = 3.5J. 

contain the same information because in this case 

|Re(0K(iK(0)|0>| « |<VK(*K(0M|. (7) 

In the tunnelling limit, where |(0|n z |l)| ~ 1, one obtains 
(0\n z (t) n z (0)\0) ~ |(0K|l)| 2 e lZl *. The correlation func- 
tions obtained from ED indeed show coherent oscillations 
with periods 2n/A — 24.5/J near the plateau centre [Fig. 
3(a)] and 27r/zi = 132. 0/J near the level crossing, where 
in addition a strong component of a higher harmonic is 
clearly evident [Fig. 3(b)]. The solid and dashed curves in 
Fig. 3, representing respectively the left- and right-hand 
sides of the two-level approximation (R), do not coincide 



because of additional components present in the correla- 
tion function of [right-hand side of (R)]. These indicate 
that |1) has significant matrix elements of n z with states 
other than |0) (Fig. 2), as a result of which the mapping 
of Na:Fe6 onto a two-level system is marginal (below). 

In stark contrast to M zz , the total-spin correlation 
function S zz (Fig. 4) shows oscillations only at the much 
higher frequency h x /2if, as expected from the symmetry 
considerations presented above. While these field-driven 
oscillations are coherent, they do not correspond to a tun- 
nelling process, where the levels involved lie below the 
height of the anisotropy barrier, as is the case in the low- 
est manifold. 



3.1.2 Spectra and Matrix Elements 

The time-dependent correlation functions shown above 
may be understood directly from the matrix elements be- 
tween the lowest-lying energy levels for the components of 
the total-spin and Neel-vector operators. The symbols in 
Fig. 5 show the numerical results for energy- level splittings 
and matrix elements of n z in the low-energy manifold of 
NaiFeg. The energy separations Z\, [Fig. 5(a)] confirm that 
an appreciable separation remains between the lowest pair 
of states and the next higher level for all fields. That this 
criterion alone is not sufficient to assess the quality of a 
two-level description is shown by the matrix elements in 
Fig. 5(b). While (l|n z |0) is indeed large for all fields, the 
matrix element (2|n z |0) is also significant at fields close to 
the level crossings. In fact at the plateau centers (2|n z |0) 
vanishes identically, and in this field regime (l|n z |0) is 
considerably larger than (?|ra z |0) for all i > 2. Our ED 
calculations confirm that the semiclassical picture of a 
tunnelling, or coherent oscillation at sub-barrier energies 
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Fig. 6. (a) Evolution with magnetic field of energy spectra for 
Na:Fee up to Ai = 5.5J. The solid lines represent ED data 
in the momentum sectors p — and p = ir (in both cases, 
only states with positive site parity appear in the energy range 
of the figure), while the dashed lines correspond to momenta 
p — 7r/3 and p = 2n/3 which are absent in the RR model. Note 
the almost linearly Zeeman-split state at Ai m h x . (b) Matrix 
element |(i|S z |0)| as a function of magnetic field, represented 
by radius of open circles. The dominant matrix element corre- 
sponds to the level with Ai ~ h x for h x > 0.7 J, while for the 
same field range the next-largest elements correspond to still 
higher excited levels. 



of the Neel vector between directions +e z to — e x , is in- 
deed appropriate here. By contrast, at the level-crossing 
fields more of the higher matrix elements of n z are ap- 
preciable [Fig. 2(b)], and in particular (2|n z |0) has 66.0% 
of the magnitude of (l|n z |0). In the semiclassical descrip- 
tion this corresponds to the Neel vector being rather less 
well localised along ±e z than for fields at the plateau cen- 
ters. The effects of the higher matrix elements with |0) 
[Fig. 5(b)] are clear in the difference between the solid 
lines in Figs. 3(a) and (b), while those of the elements 
with |1) are visible in the differences between solid and 
dashed lines in Fig. 3. 

Fig. 6(a) shows the spectrum of Na:Fe6 as a function 
of field, expanded to splittings Ai = 5.5J, which illus- 
trates both the predominance of the anticrossing between 
the second and third levels and the presence of a linearly 
evolving Zeeman-split level at Ai w h x . The total-spin ma- 
trix elements in Fig. 6(b) show the two primary features 
expected on symmetry grounds: the matrix element in the 
low-energy sector vanishes, (1|5 Z |0) = 0, and the domi- 
nant matrix element appears at the Zeeman splitting. This 
situation remains very close to the k z = limit, where the 
field-driven Zeeman transition is the only process with a 
non- vanishing matrix element, |(«|iS' z |0)| = y S{h x )/2, in 
which S{h x ) denotes the spin of the ground state at field 
h x . 



To clarify the issue of the importance of choosing n z 
and S z , we comment that we have also studied correlation 
functions of the other components of n and S. These quan- 
tities confirm that oscillations generically similar to those 
observed in M zz and S zz remain. However, the correlation 
functions of the vectors tend to show less clearly defined 
oscillations as compared to those of the z-components, 
because the transverse components couple rather more 
strongly to further levels. 



3.1.3 Rigid-Rotor Model 

Further insight may be gained into the nature of our exact 
results by comparison with those from semiclassical ap- 
proaches [jll],[l2],[2(]], which specify the conditions for the 
two-level system to provide an appropriate description of 
the low-energy spectrum of Eq. (|). Under the assump- 
tions that the spins in the ferric wheel have alternating 
(Neel) alignment and that magnon excitations may be ne- 
glected, H can be mapped to the Hamiltonian of a rigid 
rotor (RR) ||, 



N 



h • L - Nk 



2 2 

■s n y , 



(8) 



where n and L are respectively the position and angular 
momentum of a particle confined to the unit sphere. The 
operator for total spin is represented by the angular mo- 
mentum of the particle, S = L, the eigenstates of which, 
\l,m), are spherical harmonics. The term —Nk z s 2 n 2 ac- 
counts for the anisotropy potential which renders energet- 
ically favourable those spin configurations for which the 
Neel vector n is aligned with ±e z . For small k z , the eigen- 
states of Eq. (||) have almost uniform probability distri- 
bution of n in the plane perpendicular to the magnetic 
field, which corresponds to the kinetic limit of i?RR. The 
opposite limit of large k z is specified by the condition that 
the tunnel action So = Nsy // 2k z / J be very much greater 
than 1 fl^l , and it is here that a two-level description of 
quantum tunnelling of the staggered magnetisation is ap- 
propriate. 

The condition Sq > 1 may be taken to mark the on- 
set of a spin quantum tunnelling regime, in which there is 
only one pair of tunnel-split states in the low-energy sec- 
tor. In this respect a two-level approximation is marginal 
for the real materials, although should be rather better 
defined for Cs:Fes (So = 3.8) and Feio {Sq = 3.3) than for 
Na:Fe6 (So = 2.5). In this regime of intermediate k z /J the 
RR framework remains applicable, but the anisotropy en- 
ergy —Nk z s 2 n 2 , which may not be treated perturbatively, 
gives rise to significant mixing of states of differing angu- 
lar momentum /. A quantitative comparison of RR results 
with ED then requires exact diagonalisation of Eq. (||). 
This is most easily performed in the basis \l,m) in which 
L is diagonal, and the matrix elements (I' , m' \n 2 \l , m) are 
evaluated in spherical coordinates. For moderate fields, 
the unphysical states of large angular momenta {I > 15) 
may be neglected in the RR approach, and the dimension 
of the Hilbert space is then strongly reduced from that of 
the full Hamiltonian (256 compared to 46656 for Fee)- 
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Comparisons between ED of the RR model (||) and 
the exact numerical results are shown in Fig. 5. There is 
rather good general agreement in the low-energy sector, 
especially in magnitudes of Z\, and |(i|n z |0)|, but also a 
drift in field of the predicted magnetisation step positions 
from the exact result |30|. This is thought to be largely a 
consequence of neglecting magnon excitations, by which is 
meant spin misalignments within each of the sublattices, 
and can be removed by a uniform rescaling (with a fac- 
tor of 1.036 for the parameters of Na:Feg). The RR model 
is not expected to perform as well at higher energies, a 
statement which can be quantified by inspection of the 
exact spectra shown in Fig. 6(a), where energy levels cor- 
responding to the neglected momentum sectors p = 7r/3 
and 27r/3 appear at Ai > 4 J. However, the RR prediction 
of the total-spin matrix elements shown in Fig. 6(b) re- 
mains quite accurate for the energies shown, including the 
qualitative result that no matrix elements of the total-spin 
components connect |0) and |1) and the quantitative result 
for the ground-state spin S(h x ) ~ \Nh x /AJ\ appearing 
in the only large matrix element |(z|5 z |0)| s» ^JS(h x )/2 
at Ai w h x . We may summarise by remarking that, for 
<So > 1, the two-level paradigm H,plJ delivers a sim- 
ple conceptual picture of the low-energy spin dynamics 
in terms of coherent tunnelling of the staggered magneti- 
sation, and that in addition the semiclassical treatment 
based on ED of the RR approximation provides semi- 
quantitative accuracy for the physical parameters of the 
real materials. 

3.2 Cs : Fe 8 

We turn briefly to the 8-membered ferric wheel Cs:Fes. 
The effective uniaxial anisotropy which may be extracted 
from the magnetisation data for this material [H is con- 
siderably stronger than in the case of Na:Fe6 presented 
above. Numerically, dynamical simulations remain possi- 
ble for this system, despite the much larger Hilbert space. 
Fig. 7 shows the energy spectrum for a Cs:Feg system, 
again described by Eq. ([!]), with the anisotropy ratio de- 
duced IpG] from the angle-dependence of the first critical 
field [|l4[ to be k z /J = 0.0185. The situation remains qual- 
itatively similar to Na:Fe6, but has visible differences as 
a result of the larger values of N and k z /J. The stronger 
coupling between eigenstates \l,m) results in a stronger 
anticrossing of levels |1) and \2) at the plateau centers 
and a smaller maximal A. This enhanced separation of 
the two-state manifold at lowest energies, combined with 
an increased tunnel barrier, makes the semiclassical tun- 
nelling description more appropriate, as expected from the 
larger value of So. Although in Fig. 7 we show only the 
results of full diagonalisation, we have verifed that, as in 
the case of Na:Fe6, the RR model again provides quali- 
tatively similar results, albeit with a larger drift in the 
level-crossing fields. 

Fig. 8 shows the Neel vector correlation function for 
the system in two magnetic fields chosen close to a plateau 
centre and to a plateau edge. We observe rather clean, 
single-component temporal oscillations near the centre of 




Fig. 7. Evolution with magnetic field of (a) energy-level split- 
tings Ai and (b) matrix elements (i|n z |0) for the low-energy 
sector in Cs:Fes. In (b), the circles correspond to i = 1 and 
diamonds to i > 2. 
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Fig. 8. The Neel vector correlation function N Z z{t) for Cs:Fes 
in magnetic fields (a) gfiBB x = 3.4J and (b) g^sB x — 3.72J. 



the magnetisation plateau [Fig. 8(a)], where the admixture 
of higher frequencies is indeed weak. Near the level cross- 
ing [Fig. 8(b)] the second-largest matrix element, (2|n z |0), 
has 43.4% of the magnitude of (l|n z |0), demonstrating 
again the reduced validity of the two-level description at 
these fields, compared to the plateau centers. However, 
the large separation in frequency scales results in clearly 
visible contributions from both components, and thus to 
pronounced low-frequency oscillations in Af zz correspond- 
ing to coherent Neel vector tunnelling in Fig. 8(b). 



8 



A. Honecker et aL: Spin dynamics and coherent tunnelling in the molecular magnetic rings Fe6 and Fes 



4 Discussion 

4.1 Nature of mesoscopic quantum dynamics 

With the aid of the two-level analogue we obtain a con- 
ceptual picture of the microscopic nature of quantum co- 
herence. The main result of our ED study is that the cor- 
relation functions of the total spin and the Neel vector, 
S aa (t) and N aa (t) respectively, show qualitatively differ- 
ent behaviour. Although both exhibit almost single- (or 
two-) frequency oscillations for the parameters illustrated 
here, the frequencies of these oscillations are very differ- 
ent. While S zz {t) oscillates at a frequency determined by 
the (strong) magnetic field, uj = h x , JV zz (t) exhibits oscil- 
lations with a much smaller frequency u> = A. The am- 
plitude of the Neel- vector correlation function, Af zz (t) < 
1, also allows one to conclude that the oscillations in- 
deed correspond to a true quantum tunnelling. Hence, 
although most features of the thermodynamic properties 
and the ESR spectra of ferric wheels can be understood 
from a phenomenological total-spin Hamiltonian ^, R| , 
their most interesting dynamical feature, quantum coher- 
ent spin tunnelling, is not contained in such a description. 

The clearest example of spin tunnelling would be ob- 
tained by prepar ing a two- level system in the state \ip) in- 



troduced in Sec. p.l.lj For a ferric wheel this state would 
correspond to a spin configuration of the type represented 
in Fig. 1, with the Neel vector oriented in the direction 
+e z . \ip) is a superposition of energy eigenstates whose 
evolution proceeds according to the phase factors e~ lEit 
of the energy levels, without disturbance by external pro- 
cesses (decoherence), such that after every odd number 
of half cycles the degenerate state \ip) = (|0) — |l))/-\/2, 
corresponding to the Neel vector oriented along — e z , is 
achieved. 

In the microscopic analysis of the ferric wheel, we have 
demonstrated two important features in the realisation of 
such an idealised tunnelling scenario. First, it is indeed 
possible to reduce the complex physical system to an ef- 
fective two-state model, by application of a strong field in 
the ring plane. In this case, the dominant matrix element 
of the selected operator n z is that connecting the two levels 
in the low-energy manifold, such that M zz exhibits single- 
frequency oscillations at the splitting frequency A, which 
corresponds to tunnelling dynamics of the Neel vector. We 
stress in this connection that single-frequency oscillations 
depend not only on the energetic separation of the levels 
involved (from each other or from all others in the spec- 
trum), but essentially on the matrix elements. While other 
single- frequency oscillations may be found in a dissipation- 
less quantum system, for the ferric wheel only that in the 
low-energy manifold corresponds to a tunnelling process. 
The second key feature is that in this situation the difficult 
experimental step of establishing the non-eigenstate is 
not essential, given a suitable measurement of the equilib- 
rium correlation function M zz {t). We discuss in Sec. [4.3| 
below the possibilities for executing this nontrivial task. 

Finally, we comment that the ferric wheels are interme- 
diate in size between a minimal quantum system such as 
two coupled s = 1 spins with anisotropy and the classical 



limit of a large total (staggered) spin S> 1 with at least 
two degenerate minima in an energy continuum. Compar- 
ison of our exact numerical results with analytical expres- 
sions for small systems, and with the RR model, provides 
a prescription for understanding the evolution of coherent 
phenomena with system size, and the accompanying evo- 
lution of the appropriate description from microscopic to 
semiclassical. 



4.2 Decoherence 

As already discussed in Sec. [l], we have not included deco- 
herence within our numerical calculations. On the techni- 
cal level this can be effected by the inclusion of a generic 
bath coupled to the quantum system fjlf . However, with 
the aim of establishing the nature of coherent oscillations 
in the mesoscopic ferric wheel systems, we have restricted 
the present discussion to the pragmatic level on which de- 
coherence times -T" 1 significantly longer than tunnelling 
times A~ l are a prerequisite for any observation of coher- 
ent phenomena 

In the ferric wheel systems the tunnel splitting fre- 
quencies A are very large: at the centers of the magneti- 
sation plateaux, where we have shown that the clearest 
single-frequency oscillations should be present, the tunnel 
frequencies are A = 235GHz (11.4K) in Na:Fe 6 , 62GHz 
(3.0K) in Cs:Fe 8 and 45GHz (2.2K) in Fei . At low tem- 
peratures in a sufficiently pure crystal, decoherence will 
be due primarily to additional interactions not included 
in the minimal Hamiltonian of Eq. ([!]). The intrinsic de- 
coherence rate, r, for ferric wheels is most likely to be con- 
trolled by interring dipolar interactions of electron spins 
(10-50mK), and possibly by interring superexchange pro- 
cesses, whose contributions are very difficult to estimate 
but may exceed O.IK psj . Nuclear dipolar interactions 
with *H nuclei (O.lmK) and hyperfine interactions with 
57 Fe nuclei (lmK) are sufficiently weak plj that they may 
be neglected in these systems. An approximate upper limit 
for r in Na:Fe6 may be taken from the broadening of 
the lowest magnetisation step at 40mK |j32| , which was 
measured as 0.4T (0.54K), although this is considerably 
smaller than the broadenings observed at the magnetisa- 
tion steps by NMR relaxation rate measurements in Li:Fe@ 
[^3j. In any event, even with these worst-case decoherence 
estimates, the ferric wheels Feio, Cs:Fes and Na:Fe6 re- 
main by the A 3> r criterion very promising candidates 
for observation of coherent spin tunnelling. 



4.3 Detection of coherent tunnelling 

Using the above picture of mesoscopic quantum coherence 
at the microscopic level, one would like to consider those 
experiments or techniques which may be applied in order 
to detect coherent quantum tunnelling in the ferric wheels. 
Unfortunately, our prescription is at first sight rather ex- 
acting, as it requires observing time-dependent correlation 
functions of the Neel vector. An observation would then 
depend on a probe coupling to the staggered moment of 



A. Honecker et al.: Spin dynamics and coherent tunnelling in the molecular magnetic rings Fe6 and Fes 



9 



the AF rings. The closest available options for dynamical 
investigations are the local spin-raising and -lowering op- 
erations (Asf = ±1) in NMR and INS studies. We have 
argued that single-spin dynamics in the ring indeed re- 
flects the Neel vector response, and thus conclude that 
these techniques are in principle capable of revealing the 
existence of coherent tunnelling processes at frequency A. 
The matrix elements we have computed could also be used 
to fit transition intensities higher in the INS spectrum. In 
contrast to these probes, ESR is sensitive only to the total 
spin S of the ferric wheels, and thus, as shown above and 
in Ref. |2lJ], cannot detect mesoscopic Neel vector tun- 
nelling. The matrix elements for total spin indicate that 
the ESR response would be dominated by the field-driven 
Zeeman transition. 

However, technical problems arise which would make 
any observation of quantum coherence very difficult for 
the pure ferric wheels which are known. NMR measure- 
ments suffer from the very weak matrix elements coupling 
the nuclear spins to the electronic system pi], and from 
the mismatch in frequency scales between the maximal 
A and the probe frequency. Neutron scattering studies 
would require a large, deuterated single crystal. Here we 
suggest only that the most appropriate way forward would 
be to expand the range of viable experiments by consid- 
ering rings with a strongly broken symmetry, as may be 
achieved by doping with nonmagnetic (s = 0) impurities 
p5| , such as Ga in ferric wheels. These modified ferric 
wheels retain the AF properties required to satisfy the 
condition A 3> r, and in addition possess an excess, or 
tracer, spin which may be followed with a magnetic field. 
The dynamical properties of these systems are expected 
to retain the quantum tunnelling characteristics of the un- 
doped systems [p5| , with the important additional feature 
that the total spin S also reflects the coherent tunnelling 
of n, which would then become accessible by ESR We 
will address in a forthcoming publication the microscopic 
and experimental aspects of dynamics in modified ferric 
wheels. 



5 Summary 

In conclusion, we have presented a numerical and semiclas- 
sical analysis of low-energy spin dynamics and quantum 
coherent tunnelling phenomena in the molecular magnetic 
ring systems Na:Fe6 and Cs:Fes. The energy-level spectra 
and the matrix elements for total-spin and Neel-vector 
operators computed at different magnetic field values es- 
tablish the presence of quantum coherent oscillations in 
their correlation functions. Oscillations corresponding to 
coherent tunnelling of the staggered magnetisation are 
present for systems with the physical parameters of the 
ferric wheels, with the cleanest single-frequency, or two- 
level, oscillations found in the response of the Neel vector 
at applied fields in the centers of intermediate magnetisa- 
tion plateaux. These results show that, despite the small 
system size (N = 6 or 8) and small spin quantum number 
s = 5/2, the semiclassical picture of Neel vector tunnelling 
|fl2f provides a valid picture of the low-energy dynamical 



properties of ferric wheels. Experimental observation of 
mesoscopic quantum oscillations is possible with probes 
which flip local electronic spins, and would appear to be 
most feasible in studies of broken-symmetry systems ob- 
tainable by doping of the ring materials. 
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